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We study quantitatively angular correlations in the two-particle spectrum produced 
by an energetic probe scattering off a dense hadronic target with sizeable saturation 
momentum. To this end, two parton inclusive cross sections for arbitrary projectiles 
with small color charge density are derived in the eikonal formalism. Our results are 
the following: For large momenta of the observed particles, the perturbative limit with 
characteristic back-to-back correlation is recovered. As the trigger momenta get closer 
to the saturation scale Q s , the angular distribution broadens. When the momenta 
are significantly smaller than Q s , the azimuthal distribution is broad but still peaked 
back-to-back. However, in a narrow momentum range (0.5 -j- 1.5) Q s , we observe that 
the azimuthal correlation splits into a double peak with maxima displaced away from 
180°. We argue that it is the soft multiple scattering physics that is responsible for 
the appearance of this shift in the angle of maximal correlation. We also point out 
that when the physical size of the projectile is particularly small, the double peak 
structure persists in a significantly wider range of final state momenta. 



1. Introduction 



High-pr near-side and back-to-back particle correlations have become recently 
the focus of intensive study at the Relativistic Heavy Ion Collider (RHIC). 
Such correlations give access to the microscopic dynamics of two major phe- 
nomena searched for at RHIC in order to elucidate the properties of dense 
QCD matter [1-4]: the phenomenon of perturbative saturation [5-10] ex- 
pected to determine the initial condition of hadronic collisions at sufficiently 
high center of mass energy, and the phenomenon of jet quenching [11-14] estab- 
lished to suppress the high-pr hadronic yields due to final state interactions. 
To disentangle initial and final state effects, it is clearly important to com- 
pare data on nucleus-nucleus collisions to data on proton (deuteron)-nucleus 
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collisions in which final state effects are absent. 

Two-particle correlations may provide further tests of the dynamics under- 
lying final state jet quenching, since it has been argued that their angular 
dependence is sensitive to the density of the produced matter and its collec- 
tive flow [15,16]. Moreover, it has been pointed out that in an ideal liquid, 
the energy radiated off a hard final state parton would propagate in a 'sonic 
boom' [17]. This would lead to an azimuthal back-to-back correlation which is 
peaked away from 180°, and may thus serve as a characteristic signature of a 
medium with (almost) vanishing viscosity [18]. On the other hand, it has been 
suggested that initial state saturation could affect dramatically the two-gluon 
spectrum [19] by strongly reducing the angular correlation at high transverse 
momentum, leading practically to a disappearance of the back-to-back corre- 
lations of jets in nuclear collisions. Although the effect is suggested to be more 
dramatic for large nuclear projectiles, the trend should also be seen for proton 
(deuteron)-nucleus collisions. 

The main purpose of the present work is to study quantitatively saturation 
effects in the simplest process which gives rise to two-parton correlation func- 
tions. In the target rest frame, this is the emission of one gluon from a projectile 
quark qA — > qgX with the transverse momentum of quark and gluon resolved. 
These results are presented in Section 2. 

The other goal of this work is to set up a formalism for calculating differential 
multi-parton cross sections in the eikonal approximation [20]. In particular we 
shall show how to take properly account of the final state radiation by unitarily 
transforming gluonic (and in general partonic) observables with the operator 
which dresses the partons by the cloud of Weizsacker- Williams gluons. We 
show that perturbative expansion of this cloud operator precisely generates 
the perturbative radiation in the final state. 

The paper is organized as follows: In Section 2, we discuss the physics of 
two-parton correlation functions, starting from basic formulae and numerical 
results. The subsequent sections discuss the general framework for calculating 
parton correlations in the eikonal formalism and some applications. Section 3 
presents the general discussion. This formalism is then applied to parton cor- 
relations for qA — > q(k) g(p) X in Section 4, and to qA — > q g(pi) g(p 2 ) X in 
Section 5. The main results are summarized in the Conclusions. 



2. Angular correlations in qA — > q(k) g(p) X 

We consider the simplest case of a two-particle final state. The projectile 
consists of a single quark. It scatters on a hadronic target of sizeable saturation 
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momentum and produces a gluon in the final state: q A — > q(k) g(p) X . For 
a related discussion with strong emphasis on breaking of /^-factorization, see 
Refs. [21-24]. 



A. Basic formulae and definitions 



The starting point of our study is the following expression for the doubly 
inclusive spectrum qA^> q(k) g(jp) X (for derivation see Sections 3 and 4) 

dN a. C F 1 /",,_, _,-i^_,w„_^ fx - z) • z 



/" rfz dzdxe- ik - x - !p ' (z -^ 



dydkdp 7T 2 (2vr) 4 7 (x-z) 2 z 2 
x [Q(z,x,0,z)S(z,z) + S(x,0) 

-S(x, z)S(z, 0) - S(x, z)S(z, 0)] . (2.1) 



The entire information about the target nucleus is contained in two target 
averages of products of Wilson lines, 
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S(u,u)^(j^Tr[W F \u)W F (u) 
Q(u,u,z,z) = (^Tr [w F \u)W F (u)W F \z)W F (z) 



(2.2) 
(2.3) 



To model these target averages, we use the expressions [5,21] 
S(u, u) 
<2(y,x,x,y) 



exp [— v(u — u)] , 

v(x - x) + v(y - y) - v(x - y) - v(5t - y) 



v(x) 



v(x - x) + w(y - y) - v(x - y) - w(y - x) 
x exp [-u(x - x) - u(y - y)] 
^(x - y) + v(y - x) - f (x - y) - v (x - y) 
u(x - x) + u(y - y) - u(x - y) - v(y - x) 
x exp [-u(x - y) - u(y - x)] , 

1 



x2 <gw = x 2 oio log 



x 2 A 2 



+ a 



(2.4) 



(2.5) 
(2.6) 



The function v (x) has the meaning of the target gluon field correlation function 
and is directly proportional to the gluon density in the target, hence the 
logarithmic dependence on the transverse separation. 

For numerical evaluation, we shall take A = Aqcd = 0.2 GeV. The small 
regulator a = l/(x 2 A 2 ), x c = 3 GeV -1 is chosen such that the logarithm in 
(2.6) does not turn negative and that the sensitivity on the infrared cut-off 
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is negligible for sufficiently large momentum [25]. The saturation scale Q s is 
defined implicitly as 

Q 2 = Q 2 (x 2 = 1/Q 2 ) . (2-7) 



With this definition, Q 2 = 2 GeV 2 corresponds to Q^o — 0.5 GeV 2 in Eq. 
(2.6). 

In the model (2.4)-(2.6) the properties of the single inclusive gluon radiation 
spectrum qA — > qg(p)X and its small- x evolution have been studied ex- 
tensively in recent years and are well understood [25-29]. This spectrum is 
obtained from the two-particle correlation function (2.1), integrating over the 
quark momentum k, 



dN 
dy dp 



ot s Ct 



7l z 



(27T) 



/ 



dz dze _ip - ( - s) 



z ■ z 



z 2 z 2 



X 



S 2 (z,z) + 1-S 2 (z)-S 2 (z) 



(2.8) 



In the approximation (2.4)-(2.6), one finds 
dN I f , 1 a s CV z • z 



dy dp 



-I 

7T 7 



dz dz- 



x 



(2vr) 2 
1 + e- 2v{z ~* ] - 



TV 



Z 2 Z 2 

-2v(z) _ 



ip-(z— z) 



-2u(z) 



(2.9) 



which coincides with the quasiclassical expression of [30]. 

In the limit of small transverse momenta, the relevant values of the transverse 
coordinates z, z are large. In this case, the logarithm in (2.6) becomes unim- 
portant and the Gaussian approximation Q 2 (x) = Q s is justified. Formally, 
one finds in this approximation 



dN 



dy dp 



Gaussian 



J dqe-* 



/Qs 



p 2 (p - q) 



(2.10) 



Here, p 2 T^tf ^ s ^ ne typical momentum dependence of gluons radiated off 
a high-energy quark, which received a momentum transfer q. For gluons of 
small transverse momentum p, the momentum transfer q from the target is 
accumulated according to transverse Brownian motion (q 2 ) oc Q s oc A 1 / 3 . 
For a hard gluon |p| 3> Q s , however, Eq. (2.9) shows the typical power law 
oc \ In ^ characteristic for a single hard scattering off the power-law tail of 
a hard scattering center, see (2.15) below. In this way, the ansatz (2.6) for 
the saturation momentum characterizes a medium with a physically sensible 
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interpolation between single hard and multiple soft scattering as a function 
of transverse momentum transfer. We note that this is the QCD analogue of 
QED Moliere scattering theory [31], which achieves the analogous interpo- 
lation between soft multiple and single hard target momentum transfers for 
electromagnetic scattering processes. 



B. Perturbative baseline 



As a baseline for comparison of our numerical results, we collect here the 
perturbative expressions for the single and double inclusive cross sections. 
These do not include the effects of multiple rescattering which are due to the 
existence of a large saturation scale in the target. The perturbative baseline is 
obtained by identifying the single hard interaction term in the target averages 
(2.4), (2.5). This corresponds to an expansion to first order in v (x): 

5(0,11) = 1 -v(u- u) + 0(v 2 ) , (2.11) 
Q(u, u, w, w) = 1 — [v(u — u) + t>(w — w) + t>(u — w) + v(w — u) 

-v(u- w) -v(u- w)] +0(v 2 ) . (2.12) 



Inserting this into Eq.(2.1), we find 
dN _a s C F Q 2 sfi 1 



dydkdp it 2 it p 2 (k + p) 2 k 2 



(2.13) 



To obtain (2.13), one has to keep the logarithmic dependence in the scale 
entering v (x) in (2.6). Neglecting this dependence of Q s on x amounts to 
neglecting the possibly large (albeit more rare) momentum transfer due to 
the interaction with the high momentum tail of the target fields. In this case, 
<5s(x) = Q s , one finds instead 

^^g^, (p+k) , ( , 14) 



where the quark and gluon momenta are exactly balanced. This shows that as 
long as the interaction with the target can not impart a significant kick to a 
propagating parton, large momenta in the final state can only appear due to 
the large relative momentum between the quark and the gluon in the initial 
wave function of the projectile. 

Integrating Eq.(2.13) over the quark mometum k, we obtain the perturbative 
expression for the single gluon inclusive cross section 
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dN 



dy dp 



2a s C F Q 2 1 p 2 



To determine the distribution of the total recoil momentum, we integrate (2.1) 
over the relative momentum of the quark-gluon pair in the final state. Defining 

K = k + p, q = k-p, (2.16) 



we find to leading logarithmic accuracy 



= 2 ^ In 



dydK ~ 7T 2 (K 2 ) 



K 2 

A 2 



(2.17) 



This expression exhibits the typical perturbative power-law dependence for 
the total recoil momentum K above the infrared regulator scale A. This is 
consistent with the expectation, that for large recoil K, the target behaves 
like a single hard perturbative scattering center. 

We also give here the perturbative expression for the two gluon inclusive cross 
section qA — > q g(pi) g(p2) X , with gluons produced at rapidities r) and £. 
The full expression, including the effects of the saturation is given in Eq.(5.14) 
below. Its perturbative limit is 

dN a 2 2 1 (P1 + P2) 2 , , 



drjdp&dpz 7r6 v ^p 2 p 2 ( Pl + p 2 ) 2 A 2 



C. Distribution of the recoil momentum 



The total transverse momentum K of the quark-gluon pair in the final state 
traces the distribution of recoil momentum transferred from the target. We 
calculate it by integrating (2.1) over the relative transverse momentum q. 
After performing some angular integrations, one finds 

dN a,C F r°° , r 1 /^ , r** 



dz dz d(f)-Jo(Kz) 

Jo Jo Jo z 



dyKdK it 2 

x [Q(z + z,z,0,z)S(z) +S(z) 

-S{z - z)S{z) - S(z)S(z + z)] . (2.19) 



Here, we have introduced the infrared cut-off A cut to regulate the logarithmi- 
cally infrared divergent z-integral. For the discussion in this subsection, we 
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choose A cut = A, the same infrared cut-off, which regulates the perturbative 
expression (2.17). However, as discussed below, A cut has a physical interpre- 
tation as regulator of the transverse size of the incoming projectile. We will 
discuss the implications of restricing the size of the projectile in more detail 
later. 

We have evaluated expression (2.19) numerically. In Fig. 1 we plot the ratio 
of the distribution (2.19) to the perturbative expression (2.17). For momenta 
K ^> Q s well above the saturation scale, the recoil distribution (2.19) ap- 
proaches the perturbative power law oc In and is indistinguishable from 
the perturbative expression. For low momenta K < Q s , the saturated dis- 
tribution is suppressed relative to the perturbative one. This is the manifes- 
tation of the fact that a parton propagating through the saturated target is 
very unlikely to get a kick of momentum less than of order of Q s , whereas in 
the interaction with the perturbative target on the contrary low momentum 
transfer processes are very important. Finally, Fig. 1 exhibits a fairly wide 
maximum around K — (2 3) Q s , which is due to the transfer of momentum 
of order Q s to one of the propagating partons (q or g) . The general features of 
Fig. 1 are very similar to those of the nuclear modification factor calculated 
within the same model in [25,28,29]. 

( dN/KdK ) / ( dN pert /KdK ) 

1.2 | , 1 ^_L___^ I , , , , 1 

1 / : 

0.8 - / 
0.6 - / 




1 1 1 1 1 1 1 1 1 1 1 

01 23456789 10 

K/Q s 

Fig. 1. The full two-parton (quark- gluon) correlator (2.19) normalized to its 
perturbative limit (2.17) as a function of the total pair momentum K = |k + p|. 

D. Angular Correlations - large trigger momentum 

We next study the angular dependence of the two-parton correlation function 
(2.1). Shifting z — > z + x and introducing radial coordinates, we find 
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dN a s C F 1 f [a^ [xt; 

, , ,, — : — — ■ — = — o — 77; — To x dx dz dz 

dy kdk pdp d/\<pkp n (27r) J J J 

x J d(p z d(j) 2 cos(0 2 - fe) J (VA 2 + 5 2 ) 

x [Q(z + x, x, 0, z)S(z, z + x) + S(x, 0) 

-S(x, z)S(z, 0) - S{0, z)S{z + x, 0)] . (2.20) 



Here, z • z = zz cos(0 2 — 2 ), x • z = :r,2cos(0 2 ), z ■ x = zxcos(0 2 ) and the 
notational shorthands 

A = /cx cos(A0jfc p ) + px + cos(0 2 ) — cos(0 2 ) , (2-21) 

B = kx sin(A0 fcp ) — pz sin(0 2 ) + pz sin(0 2 ) . (2.22) 



The z and z integrals in equation (2.20) are logarithmically IR divergent. Reg- 
ulating these integrals by cutting off both integrations at distances greater 
than 1/A cut , we can numerically determine the correlation of the two par- 
tons as a function of their relative azimuthal angle A(p kp and their transverse 
momentum k and p. Technically, it is sufficient to regulate either one of the 
integrations z or z, to arrive at a finite result. Physically, the cut-off 1/A cut 
limits the transverse size of the incoming wave function to z < 1/ A cut in the 
amplitude and to z < 1/A cut in the complex conjugate amplitude. Thus, we 
choose to implement the cut-off in a symmetric way. 

Fig. 2(a) shows the correlation function as a function of the azimuthal angle 
for a fixed and relatively large value of the quark momentum k = 4Q S and 
values of the gluon momentum p between Q s and 4Q S . We observe that for 
larger values of p the distribution exhibits a strong back-to-back correlation. 
As p decreases, the distribution becomes wider, and finally at p = Q s it is 
practically flat. 

The main features of Fig. 2(a) can be understood in the following simple 
picture. The incoming wave function of the projectile in the present approxi- 
mation contains only one quark and up to one gluon (see Sections 3-5). The 
total transverse momentum of the incoming quark-gluon system is zero, and 
thus they are exactly correlated back-to-back. While travelling through the 
target, each parton gets a transverse kick which in most cases is close to 
Q s . However, albeit with small probability, the impinging partons can scatter 
off the perturbative high momentum tail of the target field (Moliere scatter- 
ing). The final state configurations with large transverse momentum of both, 
quark and gluon, come mainly from two sources: either from the wave func- 
tion components of the projectile where q and g have both large momentum 
and experience soft momentum transfer, or from a low momentum component 
of the incoming wave function where the quark experiences hard scattering 
off the Moliere tail and subsequently radiates a gluon in the final state. [In 
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the present calculation only the scattered quark can radiate. The radiation 
off the gluon is higher order correction in a s , since the weight of the gluon 
component in the initial state itself is of order a s .} These events are therefore 
mostly back-to-back correlated, mirroring correlations in the initial projectile 
wave function and the back-to-back nature of the final state radiation. 



(dN/pdpkdkdAc» kp )x(1/N trigger ) 



dN / pdp kdk / const. 




(dN/pdpkdkdA^ p )x(1/N trigger ) 



dN / pdp kdk / const. 



p=4 Q s , 
p=4 Q s , 



k=4 Q s 
k=3 Q' 

k=CL 




ji/2 



n 
A <l>kp 



3tc/2 




Fig. 2. (a),(c): The azimuthal dependence of the quark-gluon correlator (2.20) 
for different values of quark (k) and gluon (p) transverse momentum. The area 
under the curves is normalized to one, as required for differential rates per trigger 
particle. (b),(d): The value of the A<£ ^-integrated double differential two-parton 
yield (2.20), const. = a ^ (2 ^ . 



On the other hand, final states with a large transverse momentum of the 
quark but a gluon with p = Q s can only arise due to Moliere scattering of the 
quark, since there are no asymmetric configurations in the initial projectile 
wave function. The direction of the hard kick the quark experiences is random 
and uncorrelated with the direction of momentum of the gluon, and thus the 
angular distribution in this regime is almost completely flat, see Fig. 2(a). For 
the intermediate values of p, the distribution interpolates between these two 
extremes. 
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The angular integrated intensity (2.20) is shown in Fig. 2(b) as a function of p 
for a trigger quark momentum k = 4Q S . It exhibits a characteristic maximum 
at p = and another less sharp maximum at p — k. The maximum at p = 
originates from the momentum distribution of the incoming projectile wave 
function, which is peaked at small values of momentum as 1/p 2 . The behavior 
near this peak is unaffected by the presence of Q s , but it is affected by the 
infrared cut-off A cut , which eliminates the low relative momentum configura- 
tions from the initial wave function. The maximum at p = k can be traced to 
the pertubative expression (2.13) which (after integration over A(pk P ) diverges 
linearly at p = k. This divergence originates from the back-to-back correla- 
tions between the quark and the gluon in the projectile wave function and 
the back-to-back final state radiation. For the saturated target however, the 
final state momentum of each parton is smeared around its initial state value 
within an interval of the width ±Q S . Thus the linearly divergent peak of (2.13) 
turns into a finite peak of width Q s in Fig. 2(b). 

Figs. 2(c,d) present analogous plots, where the momentum of the outgoing 
gluon is fixed at p — 4Q S and the quark momentum is varied. Although some 
fine details are different, the overall picture is qualitatively very similar to 
Figs. 2(a,b). 

E. Angular Correlations - trigger momenta close to Q s 

We consider now the regime where the trigger momenta are of the order of 
Q s . We concentrate on the situation when the associated momenta are equal, 
p — k. In Fig. 3, we show angular correlations when the trigger momenta are 
varied between 0.625 Q s and 1.5 Q s . These plots show a structure distinct from 
the expected back-to-back correlations. The angular correlation is not peaked 
at <f) = it, but instead at = 7r — S, where 5 is clearly greater for smaller values 
of p. The dip at <fi = it disappears for trigger momenta lower than about 0.6 Q s 
and higher than about 1.4 Q s . 

The dip at n arises due to a coherent scattering effect, namely due to soft 
multiple scattering of the initial quark-gluon components of small transverse 
size. If the size of the pair is smaller than the transverse correlation length of 
the target fields, then we expect that the pair scatters as one single object as 
it sees the same target fields. It then picks up a typical soft momentum which 
is shared equally between the quark and the gluon. [In the model defined 
by (2.6), it seems reasonable to use Q s as the typical soft scale, since the 
logarithmic correction which shifts Q s $ to Q s is due to harder target gluons. 
Also, the correlation length of the target may be estimated to be of order 
l/Qs,o-] m the initial state the quark and the gluon have momenta equal 
in magnitude and opposite in direction. In our trigger momenta kinematics, 
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where the magnitudes of the trigger momenta are also equal, the momentum 
transfer from the target must be in the direction perpendicular to the initial 
momenta. We can thus estimate the maximal correlation angle when p > Q s ,o 
by the following argument. We start with the initial momenta of the quark 
and the gluon k in , — k in . In the final state k = \t in + 5k and p = — k in + 5k 
with 1 5k | = Q S fl/2 and k in • 5k = 0. For k 3> Qs,o we find the angle between 
the momenta in the final state 



= 7T — 




(2.23) 



The magnitude of the effect seen in Fig. 3 is consistent with this rough esti- 
mate. 
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So far, we have discussed how coherent scattering of incoming parton pairs 
can lead to a shift of the backward peak away of zero. To understand why the 
dip in the two-parton correlator appears only for a narrow range of trigger 
momenta, we consider two other scattering mechanisms: i) the perturbative 
hard Moliere scattering off the large momentum tails of the target, which 
leads to a back-to-back distribution, and ii) the incoherent soft scattering for 
those components of the wave function, whose transverse size is greater than 
the correlation length ~ l/Q s ,o of the target fields. Indeed, when the trigger 
momenta are much higher than Q s , the phase space in the initial state from 
which one can get to these final momenta by the coherent scattering mech- 
anism is very small. Thus the perturbative final state radiation dominates, 
which is exactly back-to-back. The phase space for coherent scattering of in- 
coming parton pairs also disappears if k tr igger is significantly smaller than Q s , 
since the coherent scattering mechanism requires enough particles in the initial 
state with momenta k initia i < k tri g ger — f3Q s with j3 - a number of order one. 
Thus, it seems natural that the coherent scattering mechanism is dominant 
only for trigger momenta around Q s . 

To further test the mechanism suggested to underly the dip in the two-parton 
correlation functions, one can exploit the dependence of the two-parton cor- 
relation function (2.20) on the infrared cut-off A cut . In the derivation of this 
expression in Section 4, this cut-off dependence can be seen to restrict the cross 
section of the final state radiation. Thus, if one artificially raises the infrared 
cut-off so that it is close to the trigger momentum, one expects that most of 
the contribution of the back-to-back correlated, final state radiation is elimi- 
nated. One expects a dip at = n for a much larger trigger momentum, with 
maxima at a position given by the simple estimate (2.23). We have performed 
various such consistency checks to further substantiate the explanation given 
above and find that the results conform with this simple picture. 

The physical interpretation to the cut-off dependence of (2.20) suggests that 
one can use it to extract some additional information from our calculation. 
As mentioned above, and as can be seen from the derivation in Section 4, the 
cut-off on the z and z integrals restricts the size of the projectile wave function 
to 1/A cut . Thus, varying the cut-off, we can probe the physics of small size 
projectiles. For example, since the scale of 0.3 Fermi is the natural size of the 
constituent quark, one may ask what is the effect of restricting the gluon in 
the wave function to be emitted no further than 0.3 Fermi from the valence 
quark. [Alternatively, one may think of this excercise as a "poor man's" way to 
mimick the doubly inclusive spectrum of a scattering of a dipole of this size.] 
In any case, raising the cut-off to this value should give a reasonable indication 
of the upper limit on the effect that one might expect. The numerical results 
with the cut-off A cut = Q s = 3.57 A are shown in Fig. 4. In accordance with 
our expectations, the effect is more pronounced. The dip at = n now persists 
up to higher values of the trigger momentum, about 2.5 Q s . It is also much 
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deeper than for the lower IR cut-off. 

dN / kdk pdp dAij) kp / const. ; A cut =3.57A 

2 
1.6 
1.2 
0.8 
0.4 


1.2 

0.8 

0.4 


0.8 

0.6 

0.4 

0.2 


0.2 

0.15 

0.1 

0.05 



tc/2 n 3ic/2 271 

Fig. 4. Angular dependence of the quark-gluon correlator (2.20) for different 
values of the trigger momentum close to Q s , const. = as ^ F . Here, the trans- 
verse size of the incoming projectile wavefunction is restricted to a small value 

1/Acut- 

This concludes our discussion of the numerical results. We now turn to the 
derivation of the basic formulae used in this section and a more general dis- 
cussion of differential cross sections in the eikonal formalism. 



3. Differential cross sections in the eikonal approximation 

Below we follow the basic formalism of [20]. Consider an energetic hadronic 
projectile impinging on a large nuclear target. The projectile is characterized 
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by a wave function, in which the relevant degrees of freedom are the transverse 
positions and color states of the partons, 

l*<n> = E ^({«„x,})|{a,,x,}). (3.1) 

{a;,Xi} 



The color index oti can belong to the fundamental, antifundamental or adjoint 
representation of the color SU (N c ) group, corresponding to quark, antiquark 
or gluon in the wavefunction. In what follows we will consider wave functions 
with a small number of partons. 

At high energy, the propagation time through the target is short, and thus 
partons propagate independently of each other. For the same reason the trans- 
verse positions of the partons do not change during the propagation. The only 
effect of the propagation is that the wave function of each parton acquires an 
eikonal phase due to the interaction with the gluon field of the target. Thus 
the projectile emerges form the interaction region with the wave function 

\*out) = S\*in) = E Wa.^n^kftlW,^}). (3.2) 

{Qi,Xi} i 

Here S is the S-matrix, and the Ws are Wilson lines 

W(xi) = Vexp{t J dz-T a At(^z-)} (3.3) 



with A + - the gauge field in the target and T a - the generator of SU(N C ) in a 
representation corresponding to a given parton. The relative phases between 
the components of the wave function change, and the state that emerges after 
the target is no longer an eigenstate of the strong interaction Hamiltonian 
(as the incoming state is assumed to be) but rather a superposition of such 
eigenstates. It is convenient to rewrite the outgoing wave function in terms of 
the second quantized operator corresponding to the gauge rotation Eq.(3.3) 

\*out) = W\*in) (3.4) 



with 



W = exp 



i / A«(x)/(x) 



(3.5) 



Here p a (x) is the color charge density operator, integrated over the rapidities 
of the projectile, and the parameters A a are the functions of the target field 
A + defined so that 



14 



exp [iT a A a (x,)] = W{^) 



(3.6) 



Various differential cross sections are given by the expectation values of gluonic 
observables 0(a,a ! ) in the outgoing wave function at t — > oo. The evolution 
between the time of scattering and the time of measurement (t — > oo) is 
unitary, and thus does not affect inclusive observables like total or inelastic 
cross sections. But it certainly does affect the measured particle spectrum 
via the final state radiation and has to be taken into account. Therefore 
expectation values of gluonic observables 0(a,cfl) are not simply given by 
(^ out \0(a, cv)\ty ou t) with \^ ut) of (3.4) since \^ ut) has to be evolved to infi- 
nite time before the average can be taken [20]. Another way to view the final 
state radiation is to realize that the state immediately after the scattering has 
a nonvanishing overlap with the incoming state \^i n ) and similar "hadronic 
states". [By hadronic states in this context we mean the states constructed 
just like the incoming state, so that they are perturbative eigenstates of the 
QCD Hamiltonian, which in the lowest perturbative order are orthogonal to 
single gluon states. We do not imply of course that within the present ap- 
proximation one has to take into account the nonperturbative hadronization 
process.] This overlap has to be subtracted to get a result which counts only 
the produced free gluons, and not those which are contained in the incoming 
wave function or the wave function of outgoing "hadronic" states but not freed 
during the scattering process. 

Consider for example the simplest possible projectile, namely the wave func- 
tion of a quark at transverse coordinate zero and color index a, which is 
"dressed" to first order in perturbation theory. The incoming quark state con- 
tains a component coming from the splitting a — > f3b, where the gluon labeled 
by adjoint index b sits at a transverse position displaced by z from its parent 
quark. To first order, the quark wave function reads 

\a D ) = \a)+ J dz «(z) Tip \(3 ; 6(z, i, 0) 

= !«>+/ <k /OO T b aP |/3 ;6(z)>. (3.7) 



Here, the second line specifies a shorthand used below, and 

f>W = (3-8) 



denotes the perturbative Weizsacker- Williams (WW) gluon field and the index 
i labels the direction in the transverse plane. For the single inclusive gluon 
cross section, one has to calculate 

{6*\af(p) af(p)|^) ; \6V) = \V out ) - £ \a D )(a D \V 0Ut ) . (3.9) 

a 
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The subtraction (3.9) properly accounts for the evolution of the scattered 
system after it emerges from the target [20], since the only "hadronic state" 
in the final state to the first order in a s is the dressed quark. Eq.(3.9) ensures 
that the observable does not include the gluons in the wave function of this 
dressed quark. 

The above example illustrates the general statement that in calculations of 
0(a,aJ), one should not count gluons that belong to the WW cloud of any 
of the fast charged partons that constitute the projectile. Rather, a "quark" 
or a "gluon" in the final state should not be thought of as a free Fock space 
quark or gluon, but the quark (or gluon) plus its WW cloud. The change of 
basis from free to dressed partons is described by the unitary "gluon cloud" 
operator C [20] 

\a D ) = C\a) (3.10) 
with 

C = Pexp(i J rfxrfz«(z-x)[a^(z,0 + af(z,0]p?(x)) , (3.11) 

where P stands for rapidity ordering, and p|(x) denotes the total charge den- 
sity operator integrated from the rapidity of the projectile to £. The expression 
(3.11) is only valid to first order in a s , but we will not need higher orders in 
the calculations presented in this paper. The observable 

CafafC* (3.12) 

counts directly those low x gluons whose wave functions are orthogonal to 
those of dressed quarks (and dressed "valence" gluons). Since the dressed glu- 
ons are eigenstates of the propagation outside the target, they can be counted 
directly in \^ ou t) without the need to account for the additional unitary evo- 
lution to infinite time. This is also true for any other gluonic observable. Thus 
the correct expression for calculating an arbitrary gluonic observable O in the 
state t — > oo is 

(V out \CO(a,a^\V out ), (3.13) 

with the gluon cloud operator C given by Eq.(3.11). One can explicitly verify 
that this procedure reproduces Eq.(3.9) in the simple example given above. 
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4. Single gluon and quark-gluon inclusive cross sections 



We now apply the formalism of Section 3 to the calculation of the one gluon 
inclusive emission cross section for an arbitrary projectile which contains a 
small number of partons, n <C l/a s . The restriction to a small number of 
partons is necessary to treat the projectile wave function perturbatively. 

We take the wave function of the projectile in the form of a valence Fock space 
dressed by the WW field 

\*in)=C\R), (4.1) 



where \R) is an arbitrary state in the free Fock space with a small number 
of large rapidity partons. The outgoing wave function after the propagation 
through the target is 

|*o« t ) = WC\R). (4.2) 



The gluon yield per unit rapidity at transverse momentum p is given by [see 
Eq. (3.13)] 



dN 
dy dp 



(* out \Ca«\p,y)at(p,y)C^ out ) . (4.3) 



To evaluate this, we act with the gluon cloud operator C on the gluon anni- 
hilation and creation operators a and cr 

Cajtf = aj- b*(p) , (4.4) 

where bf(p) is the classical WW field associated with the charge density p 
6?(M = |rfx/,(z-x)p d (x). (4.5) 



To leading order in a s , the gluon cloud operator C commutes with the charge 
density p d . We thus find 

(a, - biip))]*^) = WW\a t - bi(p))WC\R) , (4.6) 
W(w ai - bi(wp))C\R) = WC(wa, t - bi(wp))\R) . (4.7) 

We have used here the fact that the state \R) does not contain gluons with 
rapidities corresponding to the operator a, and hence it is annihilated by a. 
The shorthand notation in (4.7) stands for 
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wa* = W A b (z)a^(z), etc. 



(4.8) 



with W - the eikonal Wilson factor in the adjoint representation. We thus 
have 



dN 



dy dp (27r) 



l — J dzdze-^ z -'\R\ [w\z)bt(z,p)-bt(z,wp) 



w(z)b*(z,p)-b*(z,wp)\ \R). 



(4.9) 



Using 



w 



(z)bf(z,p) - bf(z,wp) = W A (z) J dxf t (z - x)/(x) 

- \ dx/*(z - x)W^(x)p 6 (x) , 



(4.10) 



we find 

dN cx s Cf 1 f , , t, =1 /" , , , „/ s h/ s, (z — x) • (z — x) 
dj/dp 7T 2 (27r) 2 i 7 ^ un ;/ (z-x) 2 (z-x) 2 



(\W A \z)W A (z) + lU At (x)^(x) - W^ f (z)W^(x) - W A \^)W A {z)] ah ) T , (4.11) 



where the averages of the charge density and of the products of eikonal factors 
are taken with respect to the projectile and the target wave functions, respec- 
tively. For a translationally invariant and gauge singlet but otherwise arbitrary 
target this expression reduces to the &/r-factorized form derived in [32-34]. 



A. Quark-gluon correlation function: q A — > q(k) g(p) X 



We now turn to the simplest two-parton correlation function calculable in the 
eikonal formalism: a single quark projectile of final transverse momentum k 
which shares its recoil between a gluon of transverse momentum p and the 
target. The correlation function is 



dN 



(^ out \Caf(p,y)a1(p,y)dl(k)ds(k)C^\^ out } P . 



dy dk dp 



(4.12) 



Here, d\ is the quark creation operator at projectile rapidity, and \^ ut( a )) — 
WC\a). We average over the color charge a of the incoming quark, 

{^out\0\^ out )p = ^ a ^ out {a)\0\^ out {a)) . (4.13) 
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The quark color charge density operator is 
p a (x) = 4(x)T a V /3 (x). 



(4.14) 



This charge creates the Weizsacker- Williams field according to Eq.(4.5). The 
calculation of the previous subsection is easily repeated. The action of the 
eikonal S-matrix operator XV amounts to rotating the quark and the gluon 
creation operators by appropriate eikonal factors. This amounts to attaching 
a fundamental Wilson line W^(x) at the transverse position of the quark and 
an adjoint Wilson line W^(z) at the transverse position of the gluon. We find 



af(x,y) dj(x) C^\^ out (a)) 



= - x) 7*, W#(x)W£(z) - TL W£(x) |a> . (4.15) 



In terms of this expression we can write the quark-gluon correlator (4.12) 
explicitly. In the large- N limit, the color algebra simplifies considerably and 
one obtains for the two-particle correlation integrated over impact parameter 
b 



dN 

dydkdp ~ (2tt)^ 

x [Q(z, x, x, z)S(x, z) + 5(x, x) 

-5(x, z)5(z, x) - 5(x, z)5(z, x)] . (4.16) 



Here, we use the shorthand / x = / dx. This quark-gluon correlation func- 
tion is expressed in terms of two target averages, S(z,z) and Q(z,x, x, z) 
defined in eqs. (2.2)-(2.6). Integrated over k, it reproduces the expression for 
the single gluon inclusive emission cross section given in (4.11). This can be 
checked explicitly by inserting in (4.11) the color charge density correlator 
(p a (x)p b (x))p = 2]y£ a6 £(x)<5(x) and writing the averages over adjoint Wilson 
lines in terms of fundamental ones. 



5. The two gluon inclusive cross section: nA^ g(pi) g(p 2 ) X 



In this section we calculate the eikonal cross section for production of two 
gluons. We first give the derivation for an arbitrary perturbative projectile. 
Then we specialize to the case of a single quark projectile. 

The cross section for production of two gluons with rapidities and transverse 
momenta (77, Pi) and (£, p 2 ) is 
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dN 



drjdpid£dp 2 (27r) 



-/ 

71") 4 Jzzl 



-ipi-(z-z)-ip2-( 



u - Q )(i?|Ctwt C [ a f( z , 0< (z,0 
af (u, 77) a$(u, 77)] C^C^) . (5.1) 



Our convention is such that the rapidity 77 is closer to the rapidities of the 
valence partons in the projectile, 77 ^> £. We consider 77 and £ to be sufficiently 
close, so that no evolution effects between 77 and £ have to be taken into 
account. 

For this case, the action of the cloud operator C Eq.(3.11) on the gluon field 
operator can be restricted to the two rapidities 77 and £. We can therefore 
write it schematically as 



(5.2) 



with 



C T] = exp 
C% = exp 



6f(z)[af(z,77) + af(z,77)] 
[bt(z)+5bt(z)} [af(z,0 + «f(z,0] 



(5.3) 
(5.4) 



Here, p d (x) is the charge density operator at the valence rapidity. It deter- 
mines the classical WW field bf(z) as specified in Eq. (4.5). The color charges 
produced at rapidity 77 are measured by p^(x) and constitute an additional 
contribution 8bf(z) to the WW field which affects the gluon production at 
lower rapidities £, 



5bf{z)= / x /,(z-x)pjx) 



(5.5) 



where 



p{(x) = a? (x, 77)T b d X(x, 77) , Tl = -if 



abc 



(5.6) 



To construct a state with up to two gluons, we have to expand the cloud 
operator (5.2) up to second order in p, 

C=l+7 | z {[af(z,77)+af(z,77)]6f(z) 

+ [a?(z,0 + a?W)] 
[af(z,0 + af(z,0] [6?(z) + flftz)] 



x[af(u,77) + af t (u,77)]6f(u). 



(5.7) 
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It is now a straightforward albeit tedious matter to calculate C^WC and to 
act with it on \R). [Alternatively, one can construct |\P in ) and \^ out ) from the 
cloud operator and evaluate the gluon number operator in this |\P ut) state. 
This is done in Appendix A.] Remembering that W'af(z)W = w ab (z)cii(z) 
and p a (x)W = w ab (x)p b (x), we get after some algebra 

WU«(z,Oa%u, V ) &WC\R) 

= / - xi)/j(u - x 2 ) {(w(xi) - w{z)) p(xi)} a {w(u)p(x 2 )} 6 

- /_ U{z - Xl )/,(u - x 2 ) {^(x 2 )p(x 2 )} 6 {( W ( Xl ) - w(z)) p( Xl )} a |«> 
+ ^ / i (z-u)/ j (u-xO{Mz)-^z))^(u)T 6 ^u)p(x 1 )} a | J R>. (5.8) 



With this state, the expectation value of the observable (5.1) can be calculated 
directly. It contains terms which are quadratic, cubic and quartic in the density 
of the projectile, 



dN 



7T I J zzuu 



/ zzuu 

where 

E 2 = / /(u - z) • /(z - u) /(u - x) ■ /(x - u) (p a (x)p b (x)) P 

J xx 

x({w\u)T c w(u) (w\z)-w\u)) 

(w(z) - w(u)) w t (u)T c w(u)} ab ) , (5.10) 

S 3 = / /(u - z) • /(z - Xx) /( Xl - u) • /(u - x 2 ) 

JxiXlX 2 

x ;(p c (x 1 )p e (x 2 ) P d (x 1 ))p 

x({[ W t(x 1 ) - «;t(z)][«;(z) - w(u)]T c w^u)w(x 2 )} d£ ) T 
-(p c (x 1 )p d (x 1 )p e (x 2 )) P 

x ({[^(xi) - w\z)][w(z) - w(u)]TV(u)w(u)} de ) T 

+ / /> - u) • /( Xl - z) /(u - Xx) • /(x 2 - u) 

JXlX2Xl 

x ;(p d (x 1 )p c (x 2 ) P e (x 1 ))p 

x({ W t( X2 ) u; ( ll )T e [^(z) - ^(^[^(x!) -w(z)]} C V 
-(p c (x 2 )p d (x 1 )p e (x 1 ))p 
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X 



({w\u)w(u)T e [w^(z) - w t (u)][w(x 1 ) - u?(z)]} C V 



, (5.11) 



and 



xix 2 xix 2 

x 



/(z - xi) • /(xi - z) /(u - x 2 ) • /(x 2 - u) 

( P nx 2 )p fe (x 1 )p i (x 1 )p /l (x 2 )) P 

x({[w; t (xi) - w; f (z)][w;(xi) -w(z)]} H {w t (u)w(u)}^) T 
-(p fe (x 1 )pnx 2 )p i (x 1 )p ?i (x 2 ))p 

z)][w(xi) - w(z)]} [w t (x 2 )w(u)} S ) T 



x({[ w t( Xl ; 



<{ 



x i" 1 x. 



u; t (z)][w(x 1 ) - w(z)]} W {^(ii)™^)}^ 
+ (p fc (x 1 )p 9 (x 2 )p A (x 2 )p i (x 1 ))p 

x ({[w f (xi) - w; f (z)][w;(xi) - w(z)}} M {w\x 2 )w(x 2 )Y ) T 



(5.12) 



This is an explicit function of the color charge density correlators in the pro- 
jectile, and the correlators of the eikonal factors in the target. The term S 2 
is the probability corresponding to the process when the valence component 
of \R) emits the gluon with rapidity rj, which subsequently splits into two 
gluons with rapidities r\ and £. The term S 4 corresponds to the probability of 
emission of both gluons £ and 77 from the valence component, and the term 
S 3 is the interference of these two amplitudes. These expressions are general 
and valid for any perturbative projectile. Also note that we have not used the 
formal 1/N C expansion to arrive at these expressions. In these aspects, our 
expression is more general than the one given in [21] for a dipole as projectile. 



A. Two-gluon correlations of a single quark projectile: qA^> g(pi) <?(p2) X 



We now specify to the simplest projectile - a single quark. We take the incom- 
ing quark to be at the origin of the transverse plane. The color charge density 
operator is (4. 14) and we average again over the color index of the incoming 
quark as in (4.13). The color charge density correlators are 



1 



(P a (x)p ft (y))p = -l-rt(x)<y(y) , 



(p a (x)/(y)p c (z)) P 
(p a (x)p b (y)p c (z)p d (u)) P 



2N 
1 

4iV 

2m 6 6 + 8N 



(d abe + if abe ) (d ecd + if 
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x5(x)5(y)5(z)5(u) 



(5.13) 



The two-gluon correlation function (5.1) is given explicitly in terms of these 
projectile-averaged correlators. In the large- N c limit, the expression simplifies. 
Using the £77(iV c )-identities compiled in the appendix of [20], the two-gluon 
correlation function can be expressed in terms of the target averages (2.2) and 
(2.3) of two and four fundamental Wilson lines, respectively. The final result 
takes the form 



dN 



dr)dpid£dp2 

x 



2/(z 
-/> - 



+/(*) 



(27T) 



nY L 



-ipi-(z-z)-ip2-(u-u) 



/(u) ' /(u) 



- u) • /(z - u) |Q(u,u, z,z)S'(z,z) 5(u, u) + 5 2 (u, u) 

— S(u, z)S(z, u)S(u, u) — S^u, z)S{z, u)5'(u, u)} 

- u) ■ /(z) {Q(u, u, z, z)5'(u, u)5*(z, z) + S(u)S(u)S(u, u) 

-Q(u, u, 0, z)5(z) 5(u, u) - 5(u, u)S(u, z)S(z, u) 

+Q(u, 0, z, z)5(z, z) 5(u) + S 2 (u) 

-S(u)S(z)S(z,u) - ^(u)^,^^)} 
/(z — u) {Q(u, u, z, z)S'(u, u)S'(z, z) + S'(u)S'(u)S'(u, u) 

-g(u, u, 0, z)S(z) S(u, u) - 5(u, u)5(z, u)5(z, u) 

+g(u, 0, z, z)5(z, z) S(u) + S 2 (u) 

-S(u)S(z)S(z,u) - S(u)S(z)S(z,u)} 
/(z) {g(u, u, z, z)S(z, z) 5(u, u) + S 2 (u, u) 

— Q(u, u, 0, z)S(u, u) S*(z) + ^(z, z) 

-g(u,u,z,0) 1 S(u,u) 1 S(z) + l- 1 S 2 (z)- 1 S 2 (z)}] . (5.14) 



We note that the two gluon correlation function of the perturbative qq dipole 
projectile was derived in Ref. [21] in the large N c limit. We have checked that 
in the limit of arbitrary large dipole size, the final result of Ref. [21] coincides 
with (5.14). 



6. Conclusions 



In this paper, we have shown that the framework of perturbative saturation 
provides for a mechanism which can shift the maximum strength which of 
azimuthal particle correlations away from 180°. We understand the physics of 
this phenomenon in the following way. There are two basic mechanisms by 
which a quark and a gluon can be produced in the final state. The first one 
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is the QCD analog of Moliere scattering. The quark in the projectile wave 
function scatters from the perturbative part of the target gluon field and then 
radiates a gluon in the final state. This final state radiation is predominantly 
back-to-back and generates maximal correlations at angle tt. The second mech- 
anism is the multiple scattering of the qg component of the initial state on 
the target. The final states produced by this mechanism depend strongly on 
the separation between the initial quark and gluon in the transverse plane. 
For components with large separation (small relative transverse momentum 
|k — p| <C Q s in the initial state), the quark and the gluon scatter indepen- 
dently off the target fields. This process produces uncorrelated qg pairs in the 
final state and leads to the broadening of the angular distribution. However, 
the components of the incoming state which have transverse size smaller than 
the correlation length of the target fields scatter coherently, i.e. they scatter 
effectively as a single particle. These quark-gluon components pick up a typi- 
cal soft momentum (of order Q S fl, as explained in Section 3), which is equally 
shared between the two partons. As a result, the two partons emerge from 
the interaction region with momenta shifted in the same direction by an equal 
amount of order Q s ,o/2. In the initial state, the total transverse momentum 
vanishes, and thus the momenta of the quark and the gluon are balanced. The 
equal momentum transfer therefore produces angular correlation in the final 
state which is peaked away from 180°. The maximum correlation produced 
by this mechanism for large trigger momenta k can be simply estimated to 
lie about Q s ^/\/2k away from 180°. Whether this shift in the maximal cor- 
relation angle is actually observable in the spectrum depends very much on 
the relative importance of the Moliere component and the coherent soft scat- 
tering component. Our numerical results suggest that the coherent scattering 
component is dominant for the trigger momenta around Q s . 

We also observe that the coherent scattering component, which shifts the 
maximal correlation strength away from 180°, is enhanced and dominates 
over a wider range of trigger momenta, if the projectile system has smaller 
transverse size. To establish this statement, we use the fact that our calculation 
allows us to regulate the transverse size of the incoming projectile via the cut- 
off A cut . Although it is difficult to draw firm phenomeno logical conclusions, 
we believe that there is potential for an observable effect. 

Also, it would be interesting to study how this picture discussed above is 
affected by \ow-x evolution. One may expect that since Q s increases with ra- 
pidity, the maximal correlation angle at forward rapidities should move further 
away from 180°. However, \ow-x evolution does not only affect the value of the 
saturation scale, but also changes significantly the efficiency of the target as 
a function of momentum transfer. In particular, the momentum dependence 
of the target gluon distribution decreases slower at large momenta thus effec- 
tively increasing the Moliere hard scattering component. Hence, it is not clear 
to us how the azimuthal shift of the maximal two-parton correlation strength 
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will evolve with x. 



A. Alternative calculation of two-gluon correlation function (5.14) 



In this appendix, we derive the two-gluon correlation function (5.14) following 
the formulation given in [13,20]. We start from the state in (3.1) of a 
single quark with color a and its gluon cloud, expanded in perturbation theory 
up to 0(g 2 ), 




This is the incoming state with two gluons of color b and c. The crossed line 
indicates a probability conserving virtual correction, such that the state \ip^} 
is normalized to unity up to 0(g 2 ). In this approximation, the states do not 
depend on rapidity, and the rapidity labels are suppressed in the following. 

The interaction of ) with the target results in phase shifts by eikonal Wilson 
lines, which are in the fundamental (W F ) and adjoint (W A ) representation, 
respectively. This leads to the outgoing state 



l-^/dx/(x)-/(x)) W£(0) \P) 
+t J dx/(x) (T a W F (0))^ W£(x) |#&(x)> 
-\jd*dy [/(x)/(y) (T a T d W F {Q)) a W^)Wi{y) 



/(x)/(y - x) (T e W F (0)) faw£{x) W&y] 



|/3;6(x),c(y)). (A.l) 



The quark propagates at the transverse position x g = 0, and the gluons prop- 
agate at x and y, respectively. 

According to (3.9), we have to subtract the overlap with the dressed incoming 

— * — * 

state. This leads to the state \5^f a ) which to lowest order in // consists of the 
following contributions: 
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a 



a8 



m 

Wi(y) cfyJ 

Ky 4 K<°> 
W 

p. 




Here, we use the diagrammatic shorthand of Ref. [20] in which thick lines 
denote partons propagating through the target and thus picking up eikonal 
Wilson lines in the corresponding representation. The first two graphs corre- 
spond to the case that the quark emits both gluons before the target and all 
three partons propagate through the target (thick lines). The third diagram 
can be viewed as emission after interaction since the outgoing lines do not 
carry eikonal factors (thin lines) . This third diagram, as well as the remaining 
three arise from the subtraction of the overlap in (3.9). 

A convenient expression for \5^f a ) is obtained by applying identities like 
(T b W F (0)T d ) =(T b T a W F (0)) W£(0). (A.2) 



This allows us to shift the Wilson operator W F (0) such that we can write 

= "/ dxdyfo) /(y) (T a T d W F (0)) ap [w&x) (w A (0) - W A (y)) 



-W A (0) (W A (0) - W A (y)) 6(x), c(y)> 



+ J dxdy/(x)/(y-x) (T e W F (0))^ 
-W A (x)T b a c ] |/3;&(x),c(y)>. 



T e ad W A ^)W A {y) 



(A.3) 



This expression may be written in the shorthand form, introducing an operator 
C(x, y) as follows 
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rfxrfy^(xj) |/?;&(x),c(y)) 



(A.4) 



This allows us to write the two-gluon correlation in a compact way by calcu- 
lating the expectation value of the two-particle number operator in this state, 
averaged over the incoming quark color index a, 



1 

N 



$^*al at (Pi> 0°(Pi> Oa f (p 2 , v)a(P2, V)\^c 



(2tt)< 



-ipi-(z-z)-ip 2 -(u-u) , 



(0%(z,u) 'O b a c Jz,u)) T 



be 



(A.5) 



Working out the target averages as described in [20], we obtain after some 
color algebra the result given in (5.14). 
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